Universal Chaotic Scattering on Quantum Graphs 
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We calculate the S-matrix correlation function for chaotic scattering on quantum graphs and show 
that it agrees with that of random-matrix theory (RMT). We also calculate all higher S-matrix 
correlation functions in the Ericson regime. These, too, agree with RMT results as far as the latter 
are known. We conclude that our results give a universal description of chaotic scattering. 
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Purpose. Closed quantum systems that are chaotic 
in the classical limit possess universal spectral fluctua- 
tion properties. Depending on symmetry, these coincide 
with the fluctuation properties of one of Dyson's three 
canonical random-matrix ensembles These state- 

ments, originally formulated in the form of a conjec- 
ture 0, have since been demonstrated [1, |3] with the 
help of Gutzwiller's semiclassical periodic-orbit theory 
for the two-point level correlator of general chaotic sys- 
tems and, in the case of orthogonal symmetry, for chaotic 
quantum graphs [Ij. 

For open chaotic quantum systems, the fluctuation 
properties of the scattering matrix (5'-matrix) are at is- 
sue, quantified in terms of the totality of S'-matrix cor- 
relation functions. As for closed systems, it would be 
desirable to establish (at least) the complete equivalence 
between the S'-matrix two-point correlation function for 
chaotic scattering with that of random-matrix theory 
(RMT) (calculated in Ref. [6] for the case of orthogonal 
symmetry and in Ref. for the case of unitary symme- 
try). However, we are not aware of any analytical results 
for S'-matrix correlations for chaotic scattering. 

In the present Letter, we fill that gap. For the case 
of orthogonal symmetry, we calculate the S-matrix two- 
point correlation function analytically for open chaotic 
quantum graphs. We show that our result coincides with 
that of random-matrix theory (RMT). We also cal- 
culate all higher S'-matrix correlation functions in the 
Ericson regime. These coincide with RMT results as far 
as the latter are known @, M. In that regime, the 
S'-matrix elements are 

supposed [iSIIll 

to have a Gaus- 
sian distribution. Our results show that this is the case 
only for strong absorption in all channels. We argue that 
our results are universal, i.e., apply to chaotic scattering 
in general. This establishes the complete equivalence of 
chaotic scattering and the random-matrix approach to 
scattering of Ref. . 

We focus attention on chaotic quantum graphs because 
here the semiclassical expansion is exact and scattering 
theory is particularly transparent. Chaotic scattering on 
quantum graphs was introduced in Refs. [H, [l^ where 



many of its properties were displayed with the help of 
numerical simulations, see also Refs. [3]. 

Scattering Matrix. We largely follow the developments 
of Ref. [ll]. A graph is a system of V vertices labelled 
a,l3,... that are linked by B bonds. For simplicity of 
notation we assume that every vertex a is linked by a 
single bond (a/3) to every other vertex /3 ("completely 
connected graph"). Then the number of bonds is B — 
V{V — l)/2. Our results remain valid, however, if some 
bonds are missing and B is smaller than V{V — l)/2 as 
long as the graph remains connected. We consider the 
limit V ^ 1. The lengths Lt, of all bonds b = {a/3) are 
assumed to be incommensurate but similar. That is the 
first neccessary assumption for the graph to be chaotic. 
One vertex or some or all V vertices are linked by a single 
bond each (a "lead") to infinity. The number A of such 
leads obeys 1 < A < V and defines both, the number of 
scattering channels and the dimension of the S-matrix. 

On each bond or lead, waves propagate freely with 
wave number k (the same for all bonds/leads), and 
the wave function is a linear combination of amplitudes 
expjifcx} and exp{—ikx} where x is the distance to 
one of the vertices attached to the bond/lead. For all 
bonds/leads, the coefficients of the linear combination are 
determined by boundary conditions specified in terms of 
V matrices F^") of dimension V defined for each vertex 
a. The matrix F^"^ expresses the outgoing amplitudes 
on the lead and on each of the V — 1 bonds attached to 
vertex a (written symbolically as O) in terms of the in- 
coming amplitudes on the same or any other bond or the 
lead attached to a (written symbolically as I) so that 
O = F*^")Z. With (5,^ ^ a the unitary and symmetric 
matrix F*^") has the form 



J") 



(1) 
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Here p^"'^ describes backscattering on lead a, de- 
scribes scattering from bond (a/3) to lead a or vice versa, 
and a^"^ describes scattering from bond (a/3) to bond 
(a7) or vice versa. In general, the matrix ct^"^ is symmet- 
ric and subunitary. For a vertex a without lead the first 
row and column of F*^"-' are lacking, we have F^"-' = tr^"' , 
and cr^"' is symmetric and unitary. 

Given a single incident wave in channel a only, these 
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boundary conditions completely define the total wave 
function. The amplitude of the outgoing wave in channel 
/3 is the element Sapik) of the symmetric and unitary 
scattering matrix. To write Sapik) explicitly, we define 
the symmetric subunitary block-diagonal matrix S of di- 
mension 2B — V{V — 1). Each of the V diagonal blocks 
of dimension — 1) carries one of the matrices 
a = 1, . . . , All other elements of S vanish. This de- 
fines the "vertex representation" E'^^ of S. The "bond 
representation" S^^^ is obtained by a reordering of rows 
and columns. We assign to every bond (a/3) the direc- 
tion d = -\- {d — —) \i a > P {a < respectively). We 
arrange the B bonds with positive (negative) direction in 
lexicographical order and label them consecutively with 
a running index b — 1, . . . , i?. Then every directed bond 
is uniquely defined by {b,d). The total number of di- 
rected bonds is 2B. The map (j'j^^ — > <J(ai3).(a-y) defines 
the bond representation of cr^") and, thus, the bond rep- 
resentation E(^) of E. The map T^"'' Tia/S) similarly 
defines the bond representation of the vector T. In bond 
representation we define the diagonal matrix exp{— ifc£} 
with elements Sbt'^dd' exp{— ifcL;,}. Diagonal elements in 
locations that differ only in the sign of d are pairwise 
equal. Written somewhat symbolically the iS-matrix is 

So.p{k) = 5o.0P^°'^ + {TW-'T)^p (2) 

where W — exp{— ifc£} — E^^^ Expanding in pow- 
ers of E*^^-' we obtain a simple physical interpretation of 
Eq. (21). The term containing the nth power of E*^^) is 
the sum of all semiclassical trajectories that connect the 
vertices a and /? via passage through {n + l) bonds. Each 
of the traversed bonds {hd) yields the factor exp{iA:Lfc}. 

Averages. The average over k is indicated by angular 
brackets and is taken over an interval of /c-values that 
is large compared to Tr/L^in where L^i^ is the small- 
est bond length. Because of the incommensurability of 
the Lb, that average is equivalent Q to a phase av- 
erage: For any function F we have {F[ex.p{ikLi,}]) — 
(1/(2^)) J^" dcl)bF[exp{icj,b}]. ThenEq. © implies (Sc^p) 
— 5app^°''^ [ill s-i^d can, thus, be read as S* = (5) + 
where the fluctuating part is 5*** — TW^^T- It also fol- 
lows that the average of the product of any number of 
^-matrix elements is equal to the product of the aver- 
ages (as in RMT). The RMT analogue of {8^^) = p'"' 
reads ^ (5aa) =(l-a;„)/(l -I-Xq). Here Xa depends on 
the density of resonances; such dependence is conspicu- 
ously absent for graphs where p*^"-' is independent of V . 

Swper symmetry and saddle-point approximation. The 
S'-matrix correlation function (P, Q) is defined as the av- 
erage of a product of P elements of with arguments 
k + Kp, p = 1, . . . , P and Q elements of S^* with argu- 
ments k — kg, q = 1, . . . ,Q. Without loss of generality 
we assume -P > Q > 1. Since — TW^^T, it suffices 
to work out 

( n n • (3) 



We generalize the approach of Ref. [5| (where only the 
(1,1) correlator was considered). Using supersymme- 
try @,[ll], the correlator ^ is written as the {P+Q)-fo\d 
derivative of the average of a generating function G (a su- 
perintegral) . The average over k is calculated as a phase 
average over all (j>p with the help of the colour-flavour 
transformation [16] in its most general form (for P ^ Q). 
Integrating out the original integration variables gives 

(G) = J d{Z,Z)exp{-A{Z,Z)} (4) 
where A{Z, Z) is the action 

AiZ, Z) = -STrln(l - ZZ) + isTrln(l - ZzZ^'z) 
-hisTrln(l - B+^Z-'z {bIY^Zz) . (5) 

Here STr is the supertrace. All matrices are defined in 
bond and in retarded-advanced representation. The ma- 
trix Z [Z) fills the upper (lower) non-diagonal block of 
the retarded-advanced representation, respectively, and 
is explicitly given by Zbdpts.b'd'qt's' = Stb' Zb-dptsM'qt's' , 
and correspondingly for Z. Here (6, d) are defined above, 
p (q) label the blocks needed to accomodate the P (Q) 
matrices W (W*) in the correlator (jS]), respectively, t = 
1, 2 is the time-reversal index as in Ref. 0, and s = f , b 
refers to commuting (anticommuting) variables, respec- 
tively. The matrices Z and Z"^ are related as in Ref. fs'] . 
For P > Q, Z and Z are rectangular matrices. With a = 
± distinguishing the retarded and the advanced sector, 
the matrix z is 6aa'Sbb'Sdd'Stt'Sss'iSa+Spp> exp{iKpLb/2} + 
5a~5qqi exp{jKgLf,/2}). The matrix B+ {B-) is only de- 
fined for the retarded (the advanced) sector, respectively, 

and is given by (;B+ )bdptsM'd'p't's' = f^it'^ss'^pp' (S^^j,,^, -I- 
cr|A^P'), and correspondingly for BZ^. Here cr| is the 
third Pauli matrix in superspace, and A^p'> denotes the 
source term needed to generate by differentiation of G 
the matrix element WjJ" b' d' ■ "^^^ integration measure 
in Eq. ^ is the flat Berezinian. Up to this point our 
results are exact. 

We calculate (G) using the saddle-point approxi- 
mation, putting z = 1 and A^^^ — for all j = 
p,q. As in Ref. variation of the resulting ac- 
tion with respect to Z and Z yields Z'^ = Z and, 
with E(^) = (S(^))^, the saddle-point equation (1 - 
ZZy^Z = (1 - E+ZElZ)-iE+ZEl. Here E+ is given 

by 6a+^tv^ss'^pp'^bd[h'd'^ s,nd correspondingly for E„. 
The saddle-point equation holds if ZEl — Y,\Z and if 
Y^B)-^(B)* ^ We use the vertex representation E^^) 
defined above. Eq. Q below shows that all eigenvalues 
but one of every matrix ct^"'' filling one of the diagonal 
blocks of E^^^ have magnitude unity. Except for correc- 
tions of order XjV worked out below we, therefore, have 
S(^)E(^)* = 1 and, thus, E^^^E^^)* = 1. To satisfy 
ZEl = E;Z we follow Ref. [F] and write the saddle- 
point solution Y as Sbb'Sdd'Ypts,qt's', and correspondingly 
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for Y. The matrix E^^^ is closest to the generic case if 
all eigenvalues of E*^^) differ from each other. That is the 
second neccessary assumption for the quantum graph to 
be chaotic. 

Corrections to the saddle-point action are due to devi- 
ations from z — 1, and from ~ 1. Concerning 
the former, we assume |Kp|ib ^ 1, <C 1 for all 
p,q,Lh, expand [5] z and the action A up to first order 
in Kp and kq, use E'^^'E'^^* ~ 1, and drop the source 
terms. With {cIb} = {^/T^)J2b^b the average level den- 
sity the first-order term in the exponent of Eq. 
(the "symmetry-breaking term" ) is 



SY = iTridn) STr, 



^pst 



-I- STr„ 



1 - YY 



'-qst 



1 - YY 



(6) 



The trace is only over the subspaces indicated, the matrix 
K is Sss'Stt'Spp' Hp, and correspondingly for k. Deviations 
from E^^-'E^^^* = 1 are calculated by putting z — I and 
dropping the source terms. We use Eq. ([TJ, suppress the 
index a, and take p to be real (as in RMT, that suppresses 
all elastic scattering phase shifts). Since F is unitary and 
symmetric it can be unitarily transformed into 



P 

exp{-i0i}ri/2 




exp{-i</)i}ri/2 
— pexp{— 2i(/)i} 

S^^exp{i(j)f_,} . 



(7) 

Here T = 1 — is the transmission coefficient, the indices 
/i, v run from 2 to 1^ — 1, and the phases 0i, 0^ are real 
and arbitrary. Eq. ([7]) shows that aa* differs from the 
unit matrix only in the first diagonal element which is 
1 — T. Using that for all channels a we obtain for the 
correction term in the exponent of Eq. (|4]) (the "channel- 
coupling term") 



YY 
-YY 



(8) 



Vertices not coupled to a lead do not contribute to the 
sum. In Ref. [H] it was shown that for V ^ 1 the massive 
modes do not contribute. That statement applies also in 
the present, more general case. Collecting everything we 
find 



(G) = J d{Y,Y)(^-^ exp{SY + CH} . 



(9) 



The term in big round brackets contains the source terms. 

Two-point Function. For the correlation function (1,1) 
the matrices Y and Y are both square matrices of dimen- 
sion four, and it is straightforward to work out the source 
terms in Eq. 1^. Lack of space does not permit us to 
present any details. Suffice it to say that using the trans- 
formations ti2 = Y{i - ry)-i/2, t2i = Y{1 - YY)-^/^, 
writing (SaJ) for and T^") = l-\{Saa)\^, and replac- 
ing the wave-number arguments of S by energies, the re- 
sulting expressions for {S^i^{k+K)S^^{k-K)) become for- 
mally identical to the corresponding terms in Eq. (7.23) 



of Ref. [6| for all values of the number A = 1, . . . ,V of 
channels. For the terms SY and CH in Eqs. ^ and (|S]) 
that can be checked directly. In the source terms, the 
phase cancels out. These facts establish the equiva- 
lence of the two-point functions of RMT and of chaotic 
scattering on quantum graphs. 

Ericson Regime. That regime is defined by the condi- 
tion E„T(") > 1. The cross section for chaotic scatter- 
ing is expected to display Ericson fiuctuations [l3, [HI ■ 
Numerical simulations |1^ have confirmed that expecta- 
tion. Eq. ^ allows us to determine the leading terms 
in an asymptotic expansion in inverse powers of T^"^ 
of all {P, Q) -correlation functions and, thus, the com- 
plete distribution of S'-matrix elements in the Ericson 
regime. The asymptotic terms are obtained [l3| by keep- 
ing in Eq. © only terms of lowest order in Y, Y. For 
SY + CH we obtain 

~IY1 - 2nr{d^){Kp + kq)]sTYst{YpqYqp) . 

pq ^ a ^ 

For the two-point function, the calculation pj| yields 

(A„^(/c + «)A,,(fc-«)}- ^ rM„2z7r(dR)(.c + ^) ■ 

(11) 

With r(") = 1 - and the replacement of wave 

numbers by energies, this is exactly the expression ob- 
tained for RMT in Refs. [1, [H. For the general (P, Q) 
correlation function we need to find the leading-order 
contribution to the source terms. We expand the last 
term in Eq. ([5|) with Z"^ — > F and Z Y , retaining 
only terms linear in Y and Y . (Only these are of the 
form ^pq YpqYqp which, according to Eq. (fTO| . gives the 
leading-order contribution). We need P source terms 
from B+. Expanding the exponential we keep the term 
{l/i2Pp\))[j:^^STT4{B+%YpqiBl-\Yqp)f. No two 
source terms in B^^ may have the same labels. That gives 
(l/2^)npE,STr,i((a3M(p)rp,(Stri),f,p)]. Since all 

source terms in b1_ ^ must also be different, the sum over 
q goes for P — Q over all permutations of q = 1, . . . , Q. 
Equivalently we may keep all q fixed and sum over all per- 
mutations oi p = 1, . . . ,P. For P > Q there are P — Q 
source terms A^p^ in B^_^^^ that do not have a counterpart 
in B^S^ ■ For these is replaced by E. Each of the 

resulting supertraces corresponds to one of the factors in 
the flat integration measure d{Y,Y) = Y\pqd{Ypq,Yqp). 
Therefore and because of Eq. PTI)1 . each superintegral 
factorizes into PQ terms, each factor characterized by 
the pair [p, q) of indices. In factors that do not carry 
any source terms the superintegration gives unity. The 
integration over those supertraces which carry both fac- 
tors A'-^-' and ^^'^^ yields the asymptotic form (fTTj) of the 
average of a pair of S'-matrix elements. For the P — Q 
unpaired source terms A^p^ the superintegration gives a 
non-vanishing contribution for ap = /3p only. The result- 
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ing factor is 

where again the phase (pi in Eq. ([T]) cancels out. The sum 
over q arises because in the advanced block, the matrix 
bIT^ carries the same entry E in every block labelled q. 
We suppress all arguments k for brevity. That gives 

^ p=l q=l ' sel j=l 

xEn(<]^,(-.)^'k (13) 

pcring— 1 ^ ' 

The sum with index "sel" goes over all (p^g) possibilities 
to select (P — Q) matrix elements from the first factor 
on the left-hand side. These give rise to the first product 
which vanishes unless all selected elements are diagonal. 
The remaining Q elements S^, symbolically written as 
— (kj) with j = 1, . . . , Q, appear as first factors in 

the angular brackets on the right-hand side. The sum 
with index "perm" extends over all permutations of these 
elements. Each of the terms in angular brackets on the 
right-hand side is equal to the asymptotic form dTTI) of 
the two-point function. 

Eq. (|13p gives the asymptotic form of all S'-matrix cor- 
relation functions (P, Q) and, thus, the complete distri- 
bution of the fc-dependent scattering matrix in the Eric- 
son regime. For the (2j2) correlation function it agrees 
with the result of Ref. Q. If (Saa) = for all a (strong 
absorption in all channels), the factors J- all vanish, the 
correlation functions (P, Q) vanish for P ^ Q, and for 
P = Q have the form characteristic of a Gaussian ran- 
dom process. In particular, all elements of S have a 
Gaussian distribution centered at zero, and cross-section 



fluctuations have the form predicted in Refs. [13, [HI- If 
(Saa) 7^ in some channel a, that fact and the unitarity 
constraint \Saa\ < 1 distort the Gaussian distribution of 
Saa- This is the cause of the occurrence of the factor 
J- a (Eq- ([rs]) ). The factor shows that the distortion is 
biggest for 1(5*00)1 = l/\/2- For cross-section correla- 
tion functions, interest is focussed on the (2,1) correla- 
tion function. It was first noted in Refs. J9(] that this 
function differs from zero (the result given there agrees 
aympotically with our Eq. (jl3p V Implications of that 
fact for cross-section fluctuations have been discussed in 
Ref. [H. 

Conclusions. For a quantum graph to be chaotic, the 
lengths Lf) of all bonds must be incommensurate, and 
all eigenvalues of S must differ. For chaotic scattering 
on a quantum graph, we have derived formal analytical 
expressions for all (P, Q) correlation functions of the S'- 
matrix. These were used to show that the (1,1) correla- 
tion function is identical to the one obtained from RMT, 
and to calculate for all (P, Q) explicit expressions in the 
Ericson regime. The latter agree with RMT results as 
far as these are known and yield the complete S'-matrix 
distribution function in that regime. It may perhaps be 
possible to obtain from Eq. ^ explicit expressions also 
for (2,1) and (2,2). Our results apply to every chaotic 
quantum graph since we have used wave-number aver- 
ages (not ensemble averages). 

The agreement of our results and those of RMT on 
the one hand and of the two-point function for closed 
graphs ^ and that for general closed chaotic systems [1] 
on the other, strongly suggest that our results are univer- 
sal, i.e., apply to chaotic quantum scattering in general. 
That implies that chaotic scattering and the RMT ap- 
proach to scattering are completely equivalent. 
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